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Abstract 

Jet Flavor Tomography is a powerful tool used to probe the properties of Quark Gluon Plasma formed in heavy ion collisions at 
RHIC and LHC. A new Monte Carlo model of jet quenching developed at Columbia University, CUJET, was applied to predict the 
jet flavor and centrality dependence of the nuclear modification factor R AA . The predictions for fragments / = n,D, B, e, derived 
from quenched jet flavors a-g,u,c,b in central and peripheral collisions at RHIC and LHC, exhibit novel features such as a level 
crossing pattern in R^^a^f over a broad transverse momentum range which can test jet-medium dynamics in quark gluon plasmas 
and help discriminate between current energy loss models. Furthermore, the inclusion of running coupling effects seems to change 
the jet energy dependence of the jet energy loss to a non trivial constant behavior, with a visible impact on the predictions for Raa. 



1. Introduction 

We report results of a new Monte Carlo pQCD tomographic model, CUJET Q, developed as part of the of the 
ongoing Department of Energy JET Topical Collaboration [2] effort to construct more powerful numerical codes, 
which are necessary to reduce the large theoretical and numerical systematic uncertainties that have hindered so far 
quantitative jet tomography. CUJET extends the development of the the GLV, DGLV and WHDG (3] [4j O opacity 
series approaches by including several dynamical features such as: 

(1) dynamical jet interaction potentials that can interpolate between the pure HTL dynamically screened magnetic 
and static electric screening (3j [4) limits; 

(2) the ability to calculate high order opacity corrections to interpolate numerically between N = 1 and N = oo 
analytic approximations; 

(3) full jet path proper time integration over longitudinally expanding and transverse diffuse QGP geometries; 

(4) the ability to evaluate systematic theoretical uncertainties such as sensitivity to formation and decoupling phases of 
the QGP evolution, local running coupling and screening scale variations, and other effects out of reach with analytic 
approximations; 

(5) elastic in addition to radiative fluctuating energy loss distributions; 

(6) convolution over rfs and flavor dependent pQCD invariant jet spectral density (without local in p T spectral index 
approximations); 

(7) convolution over final fragmentation, Df/ a (x, <2), as well as semileptonic decay distributions. 

The model is here applied to predict the nuclear modification factor R AA ^ a ^f(y « 0,p T , yfs,C = - 5%), for 
a variety of jet parton flavors a = g,u,c,b and final fragments f = n,D, B, e~, over a broad pr kinematic range at 
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mid-rapidity for central collisions at -\[s = 0.2,2.76 GeV as observed at RHIC |7 ] and LHC (8). Our attempt is to 
provide an improved and detailed implementation of the opacity series based theoretical models and consequently set 
more stringent constraints on the nature of the Quark Gluon Plasma and its coupling with highly energetic jets. 



2. CU JET model 

We remind that CUJET uses Monte Carlo techniques to compute finite order in opacity contributions to the jet 
medium induced gluon radiative spectrum. In the following analysis, we limit our studies to the first order in opacity; 
furthermore, we model the interaction potential between the jet and the medium with the pure HTL dynamical model. 
Neither of these approximations are shown to qualitatively alter the results. The energy loss is computed starting from 
the induced radiated gluon number per collinear light cone momentum fraction, dN/dx+, where fluctuations of the 
radiated gluon number are included via Poisson expansion (incoherent emission). Elastic energy losses, estimated 
using the Thoma- Gyulassy model, are accounted for as well, and show a contribution close to 25 - 30% of the total 
energy loss. 

The plasma is assumed inhomogeneous (Glauber profile) and nonstatic (1 + 1D Bjorken expansion), and its density 
profile is constrained solely by the initial observed rapidity density dN/dy, equal here to 1000 for ^fs = 0.2 GeV 
RHIC collisions and 2200 for -\fs = 2.76 GeV LHC collisions; we also set the initial thermalization time To = 1 fm/c 
and take T ~ 100 MeV as the characteristic temperature at which the jet decouples from the medium. 

After obtaining the relative energy loss (e) probability distribution functions P a (e', Pi, x, 0) for a parton a = g,u,c,b 
of given initial momentum p h the average over initial transverse jet production points x and directions is taken 
according to the standard binary collision density and convoluted over the pQCD initial jet flavor invariant pp cross 
sections 



dyd 2 pf dyd 2 pf \J ' \d 2 pf) dyd 1 p [ 



X,(f) 



Finally, we assume in- vacuum fragmentation and compute the observed final state nuclear modification factors. 



3. Results 

Jet quenching observables such as Raa can provide quantitative tomographic information about the QGP density 
evolution and the dynamics between the jet and the medium. In order to provide more stringent constraints on the 
energy loss models and at the same time to better discriminate among them, it is necessary to perform simultaneous 
constrained predictions of as many observables as possible. In this study, even though we will be focusing on the 
nuclear modification factor only, we will examine in detail its dependence on the center of mass energy of the collision 
(which in turn regulates the density and temperature of the plasma) and its dependence on the jet parton flavor. Other 
fundamental dependencies, such as rapidity and centrality, will be object of future studies. 

In CUJET the strong coupling constant a s is the only free parameter that we fit to the data: constraining one 
reference pr = 10 GeV point of pion R n AA - 0.2 at RHIC, we set a s = 0.3. The extrapolation to LHC is then 
parameter free, assuming that a s does not vary with yfs. The results are shown in Fig.l. The most striking feature is 
the novel inversion of the n < D < B R AA hierarchy ordering at high p T , in spite of the unmodified hierarchy of the 
parton relative energy loss (AE/E) U « (AE/E) C <?c (AE/E)^. This effect is mostly due to the steeper initial invariant jet 
distributions of c and b jets at RHIC ifTOl . The same novel feature, on a broader energy range, can be observed at LHC 
(figure on the right). It becomes then evident the importance of experimentally isolating and observing the separate 
contribution of D and B mesons, since the mass splitting between c and b jets is a particularly robust prediction of 
pQCD in a deconfined QCD medium. 

Although both the pion and non-photonic electron results now agree within the experimental error to the RHIC 
data, possibly settling the issue of the "heavy quark puzzle", the LHC extrapolation seems to systematically overpre- 
dict the quenching suffered by the jets in the plasma ifTTTl . Furthermore, more recent preliminary results published 
by the ALICE and CMS collaborations fT2l appear to indicate a steeper rise of Raa in the range of momenta of in- 
terest (30 - 100 GeV), suggesting the possibility of reduced coupling observed at LHC energies and densities fT3ll . 
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Figure 1. Illustration of jet flavor tomography level crossing pattern of nuclear modification factors versus pr at y = for n, D, B, e fragmentation 
from quenched g, u, c, b jets in Au+Au 5% at RHIC (left side) and extrapolated to Pb+Pb 5% at LHC (right side) computed with the dynamic 
CUJET1.0 model at leading N = 1 order in opacity. The opacity is constrained at RHIC, given dN/dy(RHIC) = 1000, by a fit to a reference point 
R n AuAu (PT = 10 GeV) = 0.2 setting a s = 0.3. The extrapolation to LHC assumes dN c h/drj scaling of the opacity as measured by ALICE |9). The 
D, B, e bands reflect the uncertainty due to the choice of NLO or FONLL initial production spectra. Note the possible inversion of n, D, B levels 
predicted by CUJET at high pj at LHC and a partial inversion at RHIC arising from competing dependences on the parton mass of energy loss and 
of initial pQCD spectral shapes. 



Motivated by these findings, we relaxed the effective fixed alpha approximation and utilized a one-loop order running 
coupling, parametrized as follows lfT4l : 



(* s (Qr) = <<* = 



2tt 



91n(Go/A G CD) 



(Q < Go) ; 



2tt 



9ln(Q/A QCD ) 



(Q > Go) • 



(2) 



Again we choose to keep ao as the only free parameter of the model. The choice of scale Q, of the order of 1 GeV, is 
somewhat arbitrary. To address this systematic source of uncertainty, we let it vary while fixing the parameter ao to 
fit one chosen pion R^ c (p T = 40GeV) = 0.35 point. We include running coupling effects in both the radiative and 
elastic 0311 contribution to the total energy loss. The results are shown in Fig. 2. 

Observing the figure on the left, it is evident that the overall shape of Raa across the broad range of pr under 
consideration is changed with respect to the previous fixed coupling results. Besides appreciating the more satisfactory 
agreement with data, both at LHC and RHIC (in the latter case our predictions are almost left unchanged given the 
restricted range of energies at play), it is surprising to note how the effective energy dependence itself of the energy 
loss appears to be modified (figure on the right). Assuming in fact a simplified model for the energy loss 



-e =kE 



<a-l L b p c 



(3) 



and extracting the index a(E) from our results, it seems that the pQCD InE « E l/3 - E l/4 characteristic LPM depen- 
dence of the energy loss is canceled when the running coupling effects are included. 



4. Conclusions 

The CUJET model has been applied to study the flavor and dependence of the nuclear modification factors for 
central collisions at mid-rapidity. With one free parameter (a s ) used to fit the pion data at RHIC, we have predicted 
a novel level crossing pattern of Raa for different flavors. The inclusion of running coupling effects in the model has 
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Figure 2. (Left) Pion Raa at LHC (solid lines) and RHIC (dashed lines). In black the fixed coupling CUJET results, constrained at pr = 10 GeV 
RHIC with a s = 0.3. In red the running coupling CUJET results, constrained at pj = 40 GeV LHC with ckq = 0.4. Notice the unmodified behavior 
of RHIC Raa in the range of energies 5-20 GeV, which makes the running coupling extrapolation to RHIC compatible with present data fl6l . 
Central 0% - 5% preliminary ALICE and CMS h± LHC data 1 12 | (brown and gray triangles, respectively) are compared to predictions. (Right) 
Illustration of the power law index a(E) as a function of jet energy extracted from a CUJET Bjorken expanding 5 fm brick-like simulation of light 
quark quenching (radiative only). A three point moving average is shown in the figure. 



instead improved the agreement with the recent LHC data without affecting the RHIC sector in the range of energies 
currently probed. Further data may give additional indication of running effects and lead to stronger conclusions with 
respect to the energy dependence of the energy loss itself. 
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